Role of disorder in the Mott-Hubbard transition 
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We investigate the role of disorder in the Mott-Hubbard transition based on the slave-rotor rep- 
resentation of the Hubbard model, where an electron is decomposed into a fermionic spinon for a 
spin degree of freedom and a bosonic rotor (chargon) for a charge degree of freedom. In the absence 
of disorder the Mott-Hubbard insulator is assumed to be the spin liquid Mott insulator in terms 
of gapless spinons near the Fermi surface and gapped chargons interacting via U(l) gauge fields. 
We found that the Mott-Hubbard critical point becomes unstable as soon as disorder is turned on. 
As a result, a disorder critical point appears to be identified with the spin liquid glass insulator to 
the Fermi liquid metal transition, where the spin liquid glass consists of the U(l) spin liquid and 
the chargon glass. We expect that glassy behaviors of charge fluctuations can be measured by the 
optical spectra in the insulating phase of an organic material k — (BEDT — TTF^iCuziC N)z. Fur- 
thermore, since the Mott-Anderson critical point depends on the spinon conductivity, universality 
in the critical exponents may not be found. 

PACS numbers: 71.10.-w, 71.30.+h, 71.10.Fd, 71.10.Hf 



I. INTRODUCTION 

Metal-insulator transition (MIT) is one of the most 
studied subjects in condensed matter physics. However, 
even the existence of the MIT is not convincingly proven 
at zero temperature in two spatial dimensions [(2 + 1)13], 
especially when both interaction and disorder coexist. [1] 
The Mott-Hubbard MIT has been claimed to occur at the 
critical interaction strength in the Hubbard model with- 
out disorder. [2] On the other hand, it is believed that the 
Anderson MIT does not arise in two spatial dimensions 
for the case of noninteracting electrons. [3] The common 
belief seems that the MIT in the presence of both in- 
teraction and disorder would not appear in (2 + \)D be- 
cause both interaction and disorder increase an insulating 
tendency. Recent experiments challenge this belief. [4] In 
Si metal-oxide-semiconductor field-effect transistor, and 
other high mobility semiconductor devices, unexpected 
MITs have been reported although these transitions are 
questioned to be truly quantum phase transitions owing 
to the temperature ranges in these experiments. [5] 

In the preset paper we investigate the role of disorder 
in the Mott-Hubbard MIT based on the Hubbard model. 
The main questions in this paper are (1) the nature of the 
insulating phase and (2) the nature of the MIT, where 
the nonmagnetic insulating phase is assumed to be a spin 
liquid state with a Fermi surface in the absence of dis- 
order. The core in these questions is the fate of gapped 
charge fluctuations in the Mott-Hubbard insulator when 
disorder is introduced. 

Using the slave-rotor representation of the Hubbard 
model, [6] we obtain an effective field theory for the MIT 
in the presence of disorder, given by a U(l) gauge theory 
in terms of fermionic spinons and bosonic collective ex- 
citations interacting via U(l) gauge fields. [7, 8] Disorder 
couples to charge fluctuations, and affect their dynam- 
ics severely. Using a renormalization group (RG) analy- 
sis, we argue that the Mott-Hubbard MIT turns into the 



Mott-Anderson MIT since the pure MIT critical point 
becomes unstable as soon as disorder is turned on, re- 
sulting in a new stable fixed point with a finite disorder. 
Accordingly, the resulting insulating phase with disorder 
is expected to be a Bose glass state, where gapped charge 
excitations in the Mott-Hubbard insulator are gapless in 
the presence of disorder. As a result, the U(l) spin liquid 
state would coexist with the Bose glass phase of charge 
fluctuations, thus called the spin liquid charge glass. 

The present study is expected to apply to geomet- 
rically frustrated lattices such as an organic material 
k~(BEDT~TTF) 2 Cu 2 (CN) 3 since the U(l) spin liquid 
Mott insulator is believed to appear in this material. [7] 
Our study implies that although spin dynamics is little 
affected by weak disorder, charge dynamics is severely 
modified from the Bose-Mott insulator to the Bose glass. 
This will be measured by charge spectra in optical con- 
ductivity experiments. 



II. SLAVE-ROTOR THEORY WITH DISORDER 



A. Formulation 



We consider the Hubbard model with disorder 

H = -tJ2 4c* + u EC 4c-) 2 
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where t is a hopping integral, u the strength of on-site 
Coulomb interaction, and vi a random potential intro- 
duced by disorder. 

The slave-rotor representation is utilized in order to 
treat the Hubbard u term. Because this methodology is 
well introduced in Refs. [6, 8], here we do not discuss 
the rotor representation in detail. An electron annihila- 
tion operator can be decomposed into a spin annihilation 
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operator f ia and a charge one e t6i in the following way 

c ia = 6 1 fia- (2) 

In this paper we call fi„ and e~ t6i spinon and chargon, 
respectively. Inserting this decomposition into Eq. (1), 
one can obtain 

Z = J D[f ia ,6 i ,Vi,L i ]e-S dTL , 

L = J2 ftAdr - Mi* -t'E&eW-Wfi, 

icr ija 

+ £>L? - iLi{d T 0i - ivi) + m(Li - ]T fUur)]- 

i a 

(3) 

Here ifi is a Lagrange multiplier field imposing the rotor 

constraint Li = J2a fLfic an d A* t ne chemical poten- 
tial of electrons. Physically, ip t is an effective electric 
potential associated with a charge density wave order 
parameter, and Li an electron density operator canon- 
ically conjugate to the phase field 9i, indicated by the 
term —iLid T 0i. It should be noted that Eq. (3) is just 
another representation of Eq. (I) via the transformation 
Eq. (2). Integrating over the potential field ipi and the 
density field Li in Eq. (3), and performing the gauge 
transformation Eq. (2), one can recover the Hubbard 
model Eq. (I). 

Integrating out the density variable L i7 Eq. (3) reads 

Z = J D[f ia ,e i ,ip i ]e-5 dTL , 

ia ija 

+ 4~^ ^2( d T&i -fi- ivi) 2 - 
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(4) 



Note that the random potential Vi couples to the charge 
density represented by d T 8i. Decomposing the hopping 
term by using the Hubbard-Stratonovich (HS) transfor- 
mation, one can obtain the effective Lagrangian 

L = Lq + Lf + Lg, 

L = t ^{aijffij + f3ijOt*j), 

(ij) 
ia 

~t ^ ] ifitjPijfja + fjaPij fia)> 

Le = j- ^{d T t - fi - m f 



(ij) 



(5) 



where ay and f3ij are spinon and chargon hopping order 
parameters, respectively. 



A saddle point analysis results in the self-consistent 
equations 

Cy ' fia. fia) = 1> 
a 

-up = -i(d T 0i) + 2u(^2 fUia) = 2u, 



a H = C^lfiafjc) 



Pi 



iSi „-i0, 



>• (6) 



Considering low energy fluctuations around this saddle 
point, one can set 



an = ae 



fa = fie™'* , pi = (f + a iT ,(7) 



where a — \(J2a fiafj<?)\ an d — \{e l6 i e~ l8i )\ are ampli- 
tudes of the hopping order parameters, and ay and cii T 
are spatial and time components of U(l) gauge fields. 

Inserting Eq. (7) into Eq. (5), we find an effective 
U(l) gauge theory for the Mott-Anderson transition 

L f = ^2 fia(9r - ia iT )fia 
ia 



i 

— 2ta^ cos(#j — 6i — ay), 



(8) 



where the mean field potential ip cancels the chemical 
potential at half filling in the fcrmion Lagrangian, and 
the Berry phase term Sb = i J2i Jo drd T 9i resulting from 
the mean field potential ip has no physical effects at half 
filling in the boson Lagrangian, thus safely ignored. [8] 



B. Discussion in the Mean field level 

It is interesting to note that spinon dynamics is de- 
coupled to chargon dynamics in the mean field scheme 
ignoring gauge fluctuations, governed by 

L f = E fWia - tp Yl + h - c -^ 

ia {ij)a 

L e = - w ^ 2 - 2ta J2 C0S ^ -^)-(9) 



Since the chargon Lagrangian corresponds to the quan- 
tum XY model in this level of approximation, the 
coherent-incoherent transition of the phase fields occurs 
in the absence of disorder, belonging to the XY universal- 
ity class at zero temperature. The incoherent phase with 
charge gap but no spin gap is identified with the Mott- 
Hubbard insulator of a spin liquid with a Fermi surface. 
In the spin liquid there is no coherent quasiparticle peak 
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at zero energy, and only incoherent hump is observed 
near the correlation energy ±u/2.[6] On the other hand, 
the coherent phase is understood as a Fermi liquid with 
a coherent quasiparticle peak at zero energy. [6] 

Realizing that the quantum XY model Lg in Eq. (9) is 
equivalent to the boson Hubbard model for critical phe- 
nomena, one can find its effective field theoretic action 
given by [9] 



S = J drd 2 r 
+vil' ji d T ip + w\ip\ 2 ■ 



(10) 



Here ip ~ {e lti } is the effective chargon field with its 
mass to 2 <~ u/t — (u/t) c , where (u/t) c is the critical 
strength of local interactions, associated with the Mott 
transition in the mean field level. [6] is a phenomcno- 
logically introduced parameter for local chargon interac- 
tions, v is a Gaussian random variable resulting from 
the Berry phase contribution, and w a Gaussian ran- 
dom mass originating from the random chemical poten- 
tial in the boson Hubbard model, where they satisfy 
(v(r)) = 0, (v(r)v(r')) = VS(r - r') and (w(r)) = 0, 
(w(r)w(r')) — WS(r — r'), respectively. 

To take into account the random variables, one can 
utilize the standard replica method, then obtain 



= J D^e- S , 
I •* 

1,1' J J 

^Jdrdr' Jd 2 r^ lT \ 2 \^ T ,\ 2 , (11) 



where 1,1' — 1, N are replica indices, and the limit of 
N — > is performed in the final stage of calculations. 



The pure Mott critical point (m 



V* — 0, and W* = 0) can be easily checked to be 
unstable against the presence of disorder (W ^ 0). In 
Eq. (11) the bare scaling dimension of ip is given by 
[ip] = L x l 2 owing to the quadratic derivative in the ki- 
netic energy term, where L is a length scale. Thus, we 
obtain [W] = L 2 , indicating that disorder is relevant at 
the pure Mott critical point. [V] — L° is obtained, thus 
marginal. In the approximation where is the 

flavor number of boson fields, it was shown that V is 
irrelevant, and W relevant so that the RG flow goes to 
the strong disorder regime. [10] But, recent RG calcula- 
tions exhibit a weak disorder fixed point in the case of 
V = 0, identified with the Bose glass to supcrfluid tran- 
sition instead of the Bose-Mott insulator to supcrfluid 
transition. [11] The weak disorder fixed point can be also 
shown to exist in the dual vortex formulation, where the 
quantum XY model is mapped into the scalar quantum 



electrodynamics in (2 + 1)D (QED3) in terms of vor- 
tices interacting via vortex gauge fields. A random mass 
term for the vortices is also induced by disorder, mak- 
ing it unstable the pure Mott critical point associated 
with the Bose-Mott insulator to supcrfluid transition. A 
new stable disorder fixed point appears to be identified 
with the Bose glass to superfluid transition in the dual 
formulation. [12] In this respect the spin liquid Mott in- 
sulator turns into the spin liquid Bose glass in the mean 
field level. 

Beyond the mean field level, dynamics of spinons and 
chargons is coupled via U(l) gauge fluctuations. In 
this case several important issues arise even in the ab- 
sence of disorder. One can doubt the stability of the 
spin liquid phase against U(l) gauge fluctuations, espe- 
cially instanton excitations allowed by the compactness 
of the U(l) gauge field. The present author discussed 
the confinement-dcconfincmcnt problem in the presence 
of the Fermi surface, and proposed the stability of the 
spin liquid phase when the spinon conductivity is suffi- 
ciently large. [13] In addition, the XY transition nature in 
the mean field level without disorder should be modified 
by spinon excitations. Especially, gapless spinon exci- 
tations result in dissipative dynamics of the U(l) gauge 
field. These damped gauge fluctuations are expected to 
turn the XY transition into the other. Furthermore, the 
presence of disorder makes the MIT much more complex. 



III. RENORMALIZATION GROUP ANALYSIS 

A. Boson-only effective action 

To investigate the role of spinon excitations in the 
coherent-incoherent transition of chargon fields, we ob- 
tain the effective chargon-gauge action in the continuum 
limit by integrating out spinons in Eq. (8) 



— (d T 6 -a T - iv) 2 - 2tacos(W - a) 
4it 



S e ff = J drd 2 r 
+ J d ^r^a^{q r ,i(jj n )D~l{q r ,iu) n )a v {-q r ,-iuj n ), 



(12) 



where D tlv (q r ,iu) n ) is the renormalized gauge propagator, 
given by 



D^ v (q r ,iuj n ) 



D(q r ,iui n ), 



D 1 (q r ,iuj n ) = D 1 (q r ,iuj n ) +n(q r ,iuj n ). (13) 

Here DQ 1 (q r , iuj n ) = (q 2 +tx> 2 )/g 2 is the bare gauge prop- 
agator given by the Maxwell gauge action, resulting from 
integration of high energy fluctuations of spinons and 
chargons. g is an internal gauge charge of the spinon 
and chargon. U(q r ,iuj n ) is the self-energy of the gauge 
field, given by the correlation function of spinon charge 
(number) currents. Since the current-current correlation 
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function is calculated in the noninteracting fermion en- 
semble, its structure is well known[14, 15] 



Il(q r ,iLu n ) = a(q r )\uj n \ + xQr 



(14) 



Here the spinon conductivity a(q r ) is given by cr(q r ) w 
ko/q r in the clean limit while it is a(q r ) « <jq = kgl in 
the dirty limit, where k is of order k F (Fermi momen- 
tum), and I the spinon mean free path determined by 
disorder scattering. The diamagnetic susceptibility \ i s 
given by \ ~ w^ 1 , where m/ ~ (i/?) -1 is the band mass 
of spinons. The frequency part of the kernel IL(q, ioJ n ) 
shows the dissipative propagation of the gauge field ow- 
ing to particle-hole excitations of spinons near the Fermi 
surface. 

Recently, the present author investigated the Mott- 
Hubbard MIT based on the effective chargon-gauge ac- 
tion Eq. (12) without disorder. [8] In this study we found 
that dissipative gauge fluctuations result in a new criti- 
cal point, depending on the spinon conductivity <jq that 
determines the strength of dissipation. In the limit of 
cr — > oo identified with a perfect metal of spinons, gauge 
fluctuations arc completely screened by spinon excita- 
tions, thus safely ignored. The resulting chargon action 
is nothing but the XY Lagrangian, yielding the XY tran- 
sition. On the other hand, in the limit of <jq — > consid- 
ered as an insulator of spinons, only the Maxwell gauge 
action is expected to appear from high energy contri- 
butions of spinons and chargons. The resulting chargon- 
gauge action coincides with the scalar QED 3 , yielding the 
inverted XY (IXY) transition[16] owing to gauge excita- 
tions. Varying the spinon conductivity, these two limits 
would be connected. 

Emergence of the new charged fixed point can be eas- 
ily understood from the effective gauge-only action at 
the critical point. Integrating over critical chargon fluc- 
tuations in Eq. (12), the critical gauge action can be 
obtained in a highly schematic form at the critical point 

S g = J d 2 Qr^a T (q r ,iuj n )Il(q r ,iuj n )& T (-q r , -iw n ), 

where a. T (q r ,ioj n ) represent the transverse components 
of the gauge fields. The gauge kernel H(q r ,iw n ) is given 

by 



I \ Q I 

U(q r ,iw n ) = — V?r + w n + <J o|w„|, 

where Ng is the flavor number of the chargon field, here 
Ng = 1 . The first term results from critical chargon fluc- 
tuations while the second originates from gapless spinon 
excitations near the Fermi surface. This gauge action 
can be easily checked to be scale-invariant at the tree 
level, giving the IXY fixed point in the <7q — > limit 
and the XY one in the Co — > oo limit. Thus, the spinon 
contribution characterized by the spinon conductivity <j 
connects these two fixed points smoothly. A finite con- 
ductivity causes a new critical point between the XY and 



IXY fixed points. In the present paper we examine the 
role of disorder in the new fixed point. 

Before closing this section, we summarize the effective 
boson-only action depending on the spinon conductivity 
(T in Table I. 



B. Dual vortex action with disorder 

Disorder effects produce random Berry phase to char- 
gon fields. Because the Berry phase term leads to a com- 
plex phase factor to the partition function of Eq. (12), it 
is not easy to handle the partition function in the char- 
gon representation. Duality transformation is generally 
performed to treat the Berry phase term. [17] The dual 
vortex action of Eq. (12) is obtained to be 



?„=/ 



drd 2 r 



+u(d x c)l + — (d x cf r - v(d x c) T - ia^d x c) A 
+ ~pYl J d qr^a^{q r ,i^n)D^{q r ,iuj n )a u {-q r ,-iuj n ). 

(15) 

Here $ is a vortex field, and a vortex gauge field. m v is 
a vortex mass, given by ml ~ (u/t) c — u/t with the mean 
field MIT critical point (u/t) c , and phenomcnologi- 
cally introduced parameter for local interactions between 
vortices. The random potential v plays the role of ran- 
dom magnetic fields in vortices. 

Since Eq. (15) is quadratic in gauge fluctuations o^, 
one finds the effective vortex-gauge action by performing 
the Gaussian integration for the gauge fields a M 



/ D[*, Cli ] 



-s v 



drd 2 r 



+u(d xc)?|- — (<9 x c)l - v(d x c) T 

+ J dTdT 1 d 2 rd 2 r 1 ^c t ,(r,T)K l _ lu (r - n,r- Tx)c v {n,Ti), 

(16) 

where the renormalized gauge propagator K^ir — ri, t — 
ti) is given by in energy- momentum space 

K^(q r ,i(j n ) = (Sf,,, - ^^jK(q r ,iuj n ), 



K(q ri iu n ) 



(Q 2 + ^Dlg 2 + v(q r )\un\ + XQr 

Q 2 + bJ 2 



(q*+u%)/g 2 + <T(q r )\w n \ 



(17) 



Here g is a redefined variable including the susceptibility. 
In the following we consider dirty cases characterized by 
a(q r ) = <7 . 
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TABLE I: Effective action depending the spinon conductivity 



(T() — » OO 

Spinon perfect metal 




< (To < OO 

Spinon metal 


(TO -> 
Spinon insulator 


Effective chargon action XY 
Dual vortex action QED3 


QED 3 
QED 3 


+ spinon-gauge correction 
+ spinon-gauge correction 


QED 3 
XY 



In order to take into account the random potential by 
disorder, we use the replica trick to average over disorder. 
The random magnetic field v in the vortex action Eq. 
(16) would cause 



E / / 

IM J J 



d 2 r-(d x ci) T {d x Q') 7 



for the Gaussian random potential satisfying (v(r)) = 
and (v(r)v(ri)) = 3(5(r — n) with the strength 5 of the 
random potential. However, inclusion of only this corre- 
lation term is argued to be not enough for disorder effects. 
Because the gauge-field propagator has off-diagonal com- 
ponents in replica indices, the vortex-gauge interaction 
of the order S 2 e^ generates a quartic term including the 
couplings of different replicas of vortices even if this term 
is absent initially. [12] Here e v is a vortex charge. The re- 
sulting disordered vortex action is obtained to be 



Z v = j c^]e 



-s v 



S v — / drd 2 r 
1 •* 

+u(d x cif T 



1 

Ata 



(d x c ; ) r 21 



Uv_ 

2 



+ E E I dqr c Alr , iu n ) K ^ (q r , i iv n ) c v i {- q r , 
1 «„ J Z 

-Y^J drdn J d 2 r^ lT \ 2 \^ Tl \ 2 

~E/ drdTi J d2r ^ d x ci ^ d x ci '^ ^ 



with W > 0. The correlation term induced by disorder 



l,l> 



drdri 



d 2 r^\<S>, 



has the same form with the term resulting from a random 
mass term. Eq. (18) is our starting action for studying 
the role of disorder in the Mott-Hubbard MIT. 



l/(2e 2 ), where e v is a vortex charge. In the limit of 
small anisotropy the anisotropy was shown to be irrele- 
vant at one loop level. [12] Furthermore, the correlation 
term between random magnetic fluxes is also ignored. 
In the small 3 limit this term was shown to be exactly 
marginal at one loop level. [12] To address the quantum 
critical behavior at the Mott transition, we introduce the 
scaling r — e r' and r = e r', and consider the renormal- 
ized theory at the transition point m 2 = 



S V = J2 I dT'd 2 r'[Zs,\(d^ - ie vCfll )^\ 2 
1 J 

+Z u ^\^ + ^(d'x Cl ) 2 _ 

-E/ dT ' dT ^ J d 2 r'Z w ^ lT \ 2 \^, Tl \ 2 , (19) 



where Z&, Z u , Z c , and Z\y are the renormalization fac- 
tors defined by 



$ = e * l Z£® r , c^^e ^Zicpr, 



-21 



1 e 2 



e l Z u Z (S , 2 u v 



W = e~ zl Z w Z^W r . 



(20) 



In the renormalized action Eq. (19) the subscript r im- 
plying " renormalized" is omitted for simple notation. 

Evaluating the renormalization factors at one loop 
level, the RG equations are expected to be[ll, 12, 16, 
18, 19] 



= u v + (h(a , e 2 v )e 2 v + jW^u v - pu 2 v - g(a , e 2 )e^, 
= 2W+ {h{a ,e 2 v )e 2 v -KU v ^W + r]W 2 . (21) 



dl 

du v 

~df 

dW 

~dT 



C. Renormalization group analysis 

We perform an RG analysis for Eq. (18). Anisotropy 
in the Maxwell gauge action for the vortex gauge field is 
assumed to be irrelevant, and only the isotropic Maxwell 
gauge action is considered by replacing u, 1/Ata with 



Here A, £, 7, p, k, 77 are positive numerical constants, and 
h((To, e 2 ), 3 (do, e 2 ,) are analytic and monotonically in- 
creasing functions of uo, as will be explained below. 

The first RG equation for the vortex charge can be 
understood in the following way. Integrating out critical 
vortex fluctuations, we obtain the singular contribution 
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for the effective gauge action 



QED 3 with a random mass term 



c tt (q r ,iuj n )E flv (q r ,iij n )c l/ (-q, -iu n ), 



E(q r ,iuj n ) = y/q 2 + + K(q r , iuj n ) 



Qr 



where iV„ is the flavor number of the vortex field, here 
N v — 1. The first term in the kernel S(g r ,za; n ) results 
from the screening effect of the vortex charge via vortex 
polarization, causing the — Ae 4 term in the RG equation 
while the second originates from that via spinon excita- 
tions, yielding the — (£/ao)e 4 term. The first e 2 term in 
the RG equation denotes the bare scaling dimension of 
the vortex charge in (2 + 1)D. We note that the above 
critical gauge action results in the relativistic dispersion 
lo ~ q r . 

For the second and third RG equations, unfortunately, 
we do not know the exact functional forms of h(ao,e 2 ) 
and g{(iQ 1 e 2 ,) owing to the complexity of the gauge kernel. 
Owing to the spinon contribution K(q r , iu> n ) [Eq. (17)] 
the kernel of the gauge propagator (c M ) 



D c (q r ,iuj n 



1 



Qr + 



e 2 K{q r , iu) n ) 



{Qr + W n)( e « +°0\Un\) 

should be utilized instead of the Maxwell propagator in 
calculating one loop diagrams. Note the dependence of 
the vortex charge e 2 in the effective gauge propagator. 
This gives the dependence of the vortex charge to the 
analytic functions h(a ,e 2 ) and g{a ,e 2 ). Although the 
exact functional forms are not known, the limiting values 
of these functions can be found. 



00 



In the limit of cr — > 00 the gauge kernel is reduced to 
the Maxwell propagator 



D c (q r ,iui n ) = 



{Q 2 +<) 



because gauge fluctuations are completely screened 
via spinon excitations in the perfect spinon metal, ig- 
nored and the resulting chargon action is nothing but 
the quantum XY model, causing the scalar QED3 as 
an effective vortex-gauge action. See Table I. Thus, 
M°o — ► 00, e 2 ) — > Ci and g(a — ► 00, e 2 ) — > c 2 are ob- 
tained, where C\ and c 2 are positive numerical constants. 
Then, Eqs. (21) become the RG equations of the scalar 



dp 2 

= e 2 - Ae 4 
dl v v ' 

du 



^ = (1 + Cl e 2 v + jW)u v - pu 2 v - c 2 e 4 , 

^- = (2 + c 1 e 2 v - K u v )W + r 1 W 2 . (22) 

These are formally the same as the RG equations stud- 
ied in Ref. [12], where the existence of the weak disorder 
fixed point was nicely discussed, guaranteeing the pres- 
ence of the Mott- Anderson MIT. 

In the absence of disorder (W* — 0) a stable charged 
critical point (e* 2 ^ 0) is expected to appear, associated 
with the Mott insulator to superfluid transition although 
there is a delicate issue about the existence of the charged 
fixed point when the flavor number of complex matter 
fields is one, corresponding to the superconducting tran- 
sition. This issue is well discussed in Ref. [16]. In this 
paper we assume the existence of the charged Mott crit- 
ical point. This fixed point becomes unstable as soon as 
disorder is turned on, as shown in the third RG equa- 
tion for W. A new stable fixed point is found with a 
finite disorder (W* ^ 0), identified with the Bose glass 
to superfluid critical point. [12] 



2. (TO 







In the spinon insulator of (Tq — > the chargon-gauge 
action is given by the scalar QED3 with disorder, as dis- 
cussed before. The resulting vortex action becomes the 
<I> 4 model with a random mass term since vortex gauge 
fluctuations c M are gapped owing to the presence of long 
range interactions mediated by the U(l) gauge fields a tll 
thus ignored in the low energy limit. This coincides with 
the fact that the gauge kernel D c (q r ,iuu) vanishes. As a 
result, h(a — > 0, e 2 ) — > and g(ao — » 0, e 2 ) — > are 
obtained. Accordingly, Eqs. (21) are reduced to the RG 
equations of the $ 4 theory with a random mass term 



du v 
~df 
dW 

~df 



= (l+jW)u v - pu 2 v , 
= {2-Ku v )W + r)W 2 . 



(23) 



The existence of the weak disorder fixed point can be 
shown in Eq. (23) [11, 20] although there are some pa- 
pers claiming that there is no weak disorder fixed point 
in this model. [9, 10] If only the two parameters corre- 
sponding to u v and W are considered in the RG equa- 
tions of Ref. [12], one finds that there indeed exists the 
weak disorder fixed point, describing the Mott- Anderson 
transition. One can perform the RG analysis based on 
the chargon-gauge QED 3 action instead of the vortex $ 4 
action. The chargon QED 3 with a random mass term 
is formally equivalent to the vortex QED3 with a ran- 
dom mass term if we correspond chargons, chargon gauge 



fields, and chargon random mass to vortices, vortex gauge 
fields, and vortex random mass. According to the pre- 
vious discussion in the vortex QED3, a disorder critical 
point would be found in the chargon QED 3 , implying 
that the Bose glass to superfluid transition also appears 
in this case. 



3. < (j < 00 

In the small do limit (<Jo\tu n \ « el) the gauge kernel 
is given by 



D c (q ri iu) n ) 



to 



e. ci + lu 2 



I) = c h a /el and g{<r ,e 2 v ) = 
z g are positive numerical con- 
stants. The corresponding RG equations are obtained 
to be 



thus resulting in h(ao,e: 
CgO-Q/ef,, where Ch and . 



pu 2 v 



^ = (1 + C h <T + "fW)u v 

— = (2 + c h cr - ku v )W + T]W 2 , 
al 



c„o. 



(24) 



where the RG flows of e 2 and u v , W are decoupled in 
this limit. If <r is replaced with e 2 in the a —* 00 limit, 
Eq. (24) coincides with Eq. (22). The replacement of 
<T with e 2 is justified by the fact that the above gauge 
kernel should be reduced to that in the <jq — > 00 limit. 
In this respect Eq. (24) can be considered to be a bridge 
between Eq. (22) and Eq. (23). 

Ignoring the (Tq term in the second RG equation, one 
finds the weak disorder fixed point depending on the 
spinon conductivity. [21] This fixed point coincides with 
that of Eq. (23) in the <r — > limit. As increasing ct , 
we expect that the fixed point of Eq. (24) gets close to 
that of Eq. (22) because Eq. (24) should correspond 
to Eq. (22), as discussed above. In other words, the 
Mott- Anderson critical point is expected to move from 
the disorder fixed point of the ao — > limit to that of the 
er —* 00 limit, depending on the spinon conductivity. 
This can be understood in the following way. The pure 
Mott critical points between the <To — > 00 and ao —> lim- 
its are smoothly connected by controlling the spinon con- 
ductivity, as discussed before. The presence of disorder 
makes the pure Mott critical points unstable, resulting in 
new disorder fixed points. Thus, it is natural that these 
new disorder fixed points are also connected smoothly 
through varying the spinon conductivity, as clearly shown 
in the small[21] and large ctq limits. Since the critical 
points depend on the spinon conductivity, the concept of 
universality is not applied to the Mott- Anderson transi- 
tion from the spin liquid charge glass to the Fermi liquid 
metal. 



w/t 



AI 



FL 



FL 



4 

1 



SLBG 



SLMI U/t 



FIG. 1: A schematic phase diagram in the slave-rotor repre- 
sentation of the Hubbard model with disorder 



D. Phase diagram and discussion 

We summarize our results in the schematic phase dia- 
gram Fig. 1, where SLBG is the spin liquid Bose glass, 
SLMI the spin liquid Mott insulator, FL the Fermi liq- 
uid metal, and AI the Anderson insulator. It should be 
noted that our approach cannot cover the whole range 
of the phase diagram. The regions indicated by question 
marks in Fig. 1 are beyond the scope of this theory. 
Strictly speaking, although the slave-rotor theory can 
produce meaningful physics in the Fermi liquid regime 
(u/t < (u/t) c ),[6] the RG equations in this paper would 
not be applied because chargon condensation (e l6i ) ^ 
allows only electron excitations owing to confinement 
between condensed chargons and spinons. When the 
strength of disorder becomes large, the RG equations 
would not work because the present analysis is based 
on the perturbation theory for weak disorder. Further- 
more, strong disorder decreases the spinon conductivity, 
making the spin liquid phase unstable against instanton 
excitations, as discussed before. [13] Remember that the 
spin liquid state can be stable in the sufficiently good 
spinon metal. Thus, our RG analysis can be applied to 
a limited range of the phase diagram near the MIT in 
the presence of weak disorder, marked by dotted arrow 
lines. In the clean limit (W — > 0) the Mott-Hubbard 
MIT is obtained between SLMI and FL.[6] On the other 
hand, in the small disorder limit the Mott-Hubbard MIT 
is shown to turn into the Mott-Anderson MIT between 
SLBG and FL. The chargon Mott insulator is expected 
to evolve into the chargon Bose glass as soon as disor- 
der is turned on, as discussed earlier. Since the chargon 
superfluidity appears in the presence of disorder, the re- 
sulting electronic phase may be identified with the Fermi 
liquid metal. We note that in the weak interaction limit 
u/t « (u/t) c the Anderson transition from FL to AI is 
expected to occur by increasing disorder although this is 
beyond the scope of the slave-rotor theory. 

A recent dynamical mean field theory (DMFT) study 
shows that the nonmagnetic phase with weak disorder 
is still a Mott insulator. [22] The Mott insulating phase 
in the DMFT study seems to be in contrast to our 
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claim that the paramagnetic phase is a gapless insu- 
lator of the Anderson type instead of the Mott one. 
We argue that this difference is not a contradiction be- 
cause physics of the DMFT approach differs from that of 
our approach. The paramagnetic Mott insulator in the 
DMFT study is different from the spin liquid Mott in- 
sulator in the slave-rotor theory in that (1) elementary 
spin excitations carry the spin quantum number 1 in- 
stead of the fractionalized spin 1/2, and (2) there is no 
spin-charge separation physics. As mentioned in the in- 
troduction, the spin liquid Mott insulator is possible to 
appear in the triangular lattice such as an organic ma- 
terial k - (BEDT - TTF) 2 Cu2(CN) 3 .[7} The present 
slave-rotor theory is expected to apply to the triangular 
lattice while the DMFT study would explain the square 
lattice. In this respect these two approaches see different 
systems, thus the resulting disordered insulators can be 
different. 

We should point out an important issue that the na- 
ture of the insulating phase in the boson Hubbard model 
with weak disorder is not completely understood. The 
Bose-Mott insulator was claimed at commensurate fill- 
ing instead of the Bose glass insulator. [23] One different 
thing from our vortex formulation is that the present 
vortex action includes anomalous gauge interactions re- 
sulting from the spinon contribution. 



IV. FERMION-ONLY EFFECTIVE THEORY 
WITH DISORDER 



There is an alternative way treating disorder in the 
slave-rotor representation of the Hubbard model. In the 
effective gauge Lagrangian Eq. (8) the gauge shift a,i T — > 
a iT — ivi results in 

i<7 t<7 

-tpJ2(f* a e- ia «f j(r + h.c.), 

(25) 



gauge action in the continuum limit 



(ij) 



Interestingly, the effect of disorder appears as the random 
chemical potential of spinons instead of the random Berry 
phase of chargons. Since the chargon dynamics does not 
couple to disorder directly, the glass phase is not likely 
to appear in this approach. Only the Mott insulator to 
superfluid transition is expected to occur in the chargon 
Lagrangian. 

Integrating out the gapped chargon excitations in the 
Mott insulating phase, we obtain an effective spinon- 



S f = J drd 2 r[^2(fUdr - ia T )U + ^jK^ ~ ia r)f*\ 



-vflh 



:\dxa\ : 



(26) 



where m/ ~ (i/3) -1 is a spinon band mass. The main 
question in this effective spinon action is about the role 
of random chemical potentials in the spinon dynamics. 

The role of nonmagnetic disorder in the QED3 with- 
out the Fermi surface was investigated by the present 
author. [24, 25] In contrast to the (2 + 1)D free Dirac 
theory long range gauge interactions are shown to re- 
duce the strength of disorder, and induce a delocalizcd 
state in the QED 3 . The presence of disorder destabi- 
lizes the free Dirac fixed point. The RG flow goes away 
from the fixed point, indicating localization. [26, 27] On 
the other hand, the charged fixed point in the QED 3 re- 
mains stable at least against weak randomness. A new 
unstable fixed point separating delocalized and localized 
phases is found. [24, 25] The RG flow shows that the ef- 
fects of random potentials vanish if we start from suffi- 
ciently weak disorder. The stability of the charged criti- 
cal point against weak disorder in the QED3, physically, 
results from the fact that the fermionic spinons feel the ef- 
fective dimensionality higher than two owing to the long 
range gauge interactions at the charged critical point, 
thus killing the effects of weak disorder. [25] 

In the spinon-gauge critical theory Eq. (26) [13] a 
similar result is expected. Deconfined spinons near the 
Fermi surface would remain delocalized at least against 
weak randomness owing to long range gauge interactions 
while noninteracting spinons without gauge interactions 
are localized by random potentials according to the scal- 
ing theory. [3] However, it should be considered that the 
presence of nonmagnetic disorder reduces the spinon con- 
ductivity (To- Thus, even if the charged fixed point can 
be stable against weak disorder in the case of noncom- 
pact U(l) gauge fields, the fixed point can be unstable 
against instanton excitations owing to the reduction of 
the conductivity. [13] As mentioned earlier, the spin liquid 
phase can be stable when the spinon conductivity is suffi- 
ciently large. In the spinon bad metal the spinons would 
be confined owing to the presence of disorder. This is 
the main reason why the application of the present slave- 
rotor formulation should be limited within weak disorder. 

The above discussion seems that the spin liquid Mott 
insulator remains stable against weak randomness in con- 
trast to the emergence of the spin liquid Bose glass in 
the first treatment. [28] To interpret this inconsistency 
between the two approaches in a consistent manner, we 
claim that the U(l) spin liquid of spinons is stable against 
weak randomness, but the Bose-Mott insulator of char- 
gons is not. The resulting insulator is identified with the 
spin liquid charge glass. 
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V. SUMMARY 

In the present paper we examined the role of disorder 
in the Mott- Hubbard metal-insulator transition based on 
the slave-rotor formulation of the Hubbard model. In 
this representation the Mott-Hubbard insulator is under- 
stood as the spin liquid Mott insulator in terms of gapless 
spinons and gapped chargons interacting via U(l) gauge 
fields. We found that the Mott-Hubbard critical point be- 
comes unstable as soon as disorder is turned on, resulting 
in a disorder critical point interpreted as the spin liquid 
glass insulator to the Fermi liquid metal transition. The 
glassy behaviors of charge fluctuations [9] can be mea- 
sured by the optical spectra in the insulating phase of an 



organic material k - (BEDT - TTF) 2 Cu 2 (CN) 3 . Fur- 
thermore, since the Mott- Anderson critical points depend 
on the spinon conductivity, universality in the critical ex- 
ponents may not be found. Last, we open the possibility 
that the spin liquid Mott insulator may survive against 
weak randomness. [23] 
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